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ABSTRACT 

We analysed the IGR J 16465^507 Burst Alert Teelescope survey data collected during the 
first 54 months of the Swift mission. The source is in a crowded field and it is revealed through 
an ad hoc imaging analysis at a significance level of ^14 standard deviations. The 15-50 keV 
average flux is '^ 3 x 10~^^ erg cm^^ s^^. The timing analysis reveals an orbital period 
of 30.243±0.035 days. The folded light curve shows the presence of a wide phase interval 
of minimum intensity, lasting ^ 20% of the orbital period. This could be explained with a 
full eclipse of the compact object in an extremely eccentric orbit or with the passage of the 
compact source through a lower density wind at the orbit apastron. The modest dynamical 
range observed during the BAT monitoring suggests that IGR J 16465— 4507 is a wind-fed 
system, continuously accreting from a rather homogeneous wind, and not a member of the 
Supergiant Fast X-ray Transient class. 
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High mass X-ray binaries (HMXBs), stellar systems composed of 
a compact object and an early- type massive star, are traditionally 
divided in two subclasses (e.g. Ivan Paradiislll995l , and references 
therein), depending on the nature of the high mass primary and, 
consequently, the different mass-transfer and accretion mechanism. 
On one side are the systems with main sequence Be primaries (Be- 
HMXBs). They are generally wide (Porb -^ 10 d) eccentric (eccen- 
tricity e ~ 0.3-0.5) systems in which the primaries are not fill- 
ing their Roche lobe, and accretion onto the compact object occurs 
from the equatorial region of the rapidly rotating Be star. Most of 
these systems are highly variable: in some of them recurrent out- 
bursts are observed caused by an enhanced rate when the compact 
star passes close to the Be star. On the other side are the systems 
with an evolved OB supergiant primary (sgHMXB). Their periods 
are shorter (Porb ^ 10 d) and their orbits more circular than in 
Be-HMXBs. They are powered either by a geometrically thin ac- 
cretion disc or by the strong radiation-driven stellar wind, depend- 
ing on whether the primary fills its Roche lobe or not. Their X-ray 
emission is bright and persistent. 

Recently, this rather clear-cut picture was made more struc- 
tured with the INTEGRAL observations of the Galactic plane. Two 
additional classes were added to the classical OB primary HMXBs: 
the highly absorbed persistent systems (Walter et al. 2004, 2006) 
and the supergiant fast X-ray transients ( SFXTs, ISguera et al.l 



X-rays: individual: IGR J16465-4507. 



former are characterized by orbital and spin periods consistent with 
those observed in wind-accreting systems, but a much higher ab- 
sorbing column density. The latter are transient sources showing a 
large dynamic range of 3-5 orders of magnitude with sporadic out- 
bursts (which however are significantly shorter than those of typi- 
cal Be-HMXBs), characterized by bright flares lasting up to days 

with peak luminosities of 10' 

iRomano et aI.ll2009l : [sidoU et al.ll2009l) 



-lO^'^ erg s"^ ( iSgueraetal., 200: 






I2OO5I : iNegueruela eTZIbOOd : ISmitlj|2004l : lin'tZandl2005l) . The 



The Bu rst Alert Tele scope (BAT, iBarthelmv eTa L 2005) on 
board Swift JGehrels et al.ll2004) is performing a continuous cov- 
erage of the hard X-ray sky (50 to 80% of the sky every day). 
This allowed the detection o f many of the new INTEGRAL 
HMXBs (e.g. ICusumano et alj 12009) and the collection of their 
long term light curves. In this Letter we analyse the hard X-ray 
data collected during the first 54 months of Swi/f-BAT sky mon- 
itoring using data in the region of the IGR J16 465-4507. This 
source was discovered by INTEGRAL in 2004 (Lutovin ov et al.l 
2004) and X-ray activity was observed with IBIS/ISGRI start- 
ing on September 6, at a flux level of 8.8 ± 0.9 mCrab (I8- 
60 keV), followed by a flare (up to 28 mCrab) on September 7. 
Follow-up o bservations with XMM-Newton revealed pulsations 
at 228±6 s ( Lutovin ov et al.l 120050 and allowed the identifica- 
tion of the optical counterpa rt with 2MASS J16463526-4507045 
jZurita Heras & Walteiil2004h . This was classified as a BO. 5 lb su- 
pergiant at a distance of ~ 8 kpc iNegueruela et alj 20071) or as 
a 09 . 5 la supergiant at a distance of Q.Slg^y^ kpc jNespoli et al.l 
120081) . The supergiant nature of the companion, combined with the 
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Figure 1. Top: 15-50 keV significance map in the neighboihood of 
IGR J16465— 4507. Tlie color-bar r epresents the significance lev els. The 
cyan star marks the XMM position iZurita Heras & Walter! 1 2004 of IGR 
J16465— 4507. The white circles are centered on the position derived by 
fitting the significance profile extracted along the white line with two Gaus- 
sians plus a constant; their radius corresponds to the 90% error on the po- 
sition. Bottom: Significance profile extracted along the white line in the 
above map. The plot shows the data (stars) and the best fit model (blue line, 
the sum of two Gaussian profiles plus a constant value). The higher peak 
(red line) corresponds to IGR J16479— 4514, the lower peak (green line) 
corresponds to IGR J16465-4507. 



observed hard X-ray variability dLutovinov et alj|2004h and the X- 
ray spectral distribut ion, modeled by a hard power law with pho- 
ton index 1.0±0.52 JLutovinov et al. 2005), suggested a classifi- 
catio n of this source as a S FXT (Negueruela et al. 2006). How- 
ever, IWalter & Zurita Heraj jIoOTh . based on INTEGRAL mea- 
surements, suggested that IGR J16465^507 is likely a classical 
supergiant HMXB, with an average flux just below the IBIS/ISGRI 
sensitivity undergoing sporadic long periods of enhanced activity. 
This Letter is organized as follows. Section 2 describes the 
BAT data reduction and the imaging analysis. Section 3 reports on 
the timing analysis. Sect. 4 describes the analysis of the pointed soft 
X-ray observation with 5wi/f-XRT. In Sect. 5 we briefly discuss our 
results. Errors are at 90 % confidence level, if not stated otherwise. 



2 IMAGING ANALYSIS 

We analysed the BAT survey data of the first 54 months of the 
Swift mission. We retrieved the raw data from the HEASARC 

S iublic archivgH and processed them with a dedicated software 
Segreto et alj 1201(1) that performs screening, mosaicking and 
source detection on data from coded mask instruments. The code 
also produces light curves for any given sky position. 

Figure [l] (top panel) shows the 15-50 keV significance 
sky map (exposure time of 17.7 Ms) in the direction of 
the sgHMXB IGR J 16465 - 4507. The position of the source 
jZurita Heras & Waited I2004I) is marked with a cyan cross. The 
source is in a crowded field, at 16.8 arcmin from the SFXT 
IGR J16479-4514 and 30.8 arcmin from the low mass X-ray bi- 
nary GX340-I-0 and for this reason it is not detected with our code. 
However, the significance profile extracted along the line crossing 
the position of IGR J16479-4514 and IGR J16465-4507 shows 
an asymmetric shape that can not be fitted with a constant plus a 
single Gaussian (the root mean square deviation is 48.8) and sug- 
gests the contribution of more than one source. The profile is in- 
stead well fitted with a constant plus a double Gaussian (the root 
mean square deviation is 3.2 and the data points are uniformly 
distributed around the model) peaking at the position of the two 
sources, with the weaker source (IGR J16465— 4507) at the signif- 
icance of 13.7 standard deviations (Fig.[T] bottom panel). 

Moreover, the presence of IGR J16465— 4507 is confirmed by 
the analysis of the sky map produced using data selected in time 
intervals wh en the nearby source IGR J16479— 4514 was in eclipse 
JBozzo et a l. 2008; Jain et al. 2009). In this map (3.4 Ms exposure) 
IGR J16465-4507 is detected with a signal-to-noise S/N~ 8. 



3 TIMING ANALYSIS 

We extracted the light curve of IGR J16465-4507 in the 15-50 
keV energy range with the maximum time resolution allowed by 
the data {^^ 300 s), su btracting the contam ination of the nearby 
sources as detailed in ISegreto et alj OOlO h. The photon arrival 
times were corrected to the Solar system barycenter (SSB) by using 
the task EARTH2SUN. 

We looked for evidence of orbital periodicities in the BAT 
data. A folding technique was applied to the baricentered arrival 
times by searching in the 0.5-100 d period range with a step reso- 
lution of P^ /{N AT), where A'^ = 16 is the number of phase bins 
and AT (140,213,559.0 s) is the data time span. The average rate 
in each phase bin was evaluated by weighting the light curve rates 
by the inverse square of the corresponding statistical error: 



R, 






(1) 



where Rj is the average rate in the j-th phase bin (j=l,16) of the trial 
profile, n are the rates of the light curve whose phase fall into the j- 
th phase bin and er^ are the corresponding statistical errors. The er- 
ror on Rj is (y^'E 1/er'j)^^ . This procedure, adopted to deal with 
the large span in er^, is justified by the fact that the data are back- 
ground dominated. Figure |2] (a) shows the resulting periodogram. 
We find significant evidence for periodicity (x ~ 251.8) at a pe- 
riod of Porb = 30.243 ± 0.035 days (the two features at higher pe- 
riods are multiple of Porb) where the error is the period resolution. 
We observe that the average x in 'h^ periodogram is far from the 

^ http://heasarc.gsfc.nasa.gov/cgi-binAV3Browse/w3browse.pl 
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average value expected for white noise (TV — 1) and increases with 
increasing trial periods. As a consequence, the x^ statistics cannot 
be applied to evaluate the significance of the detected periodicity. 
In order to have an estimate of the significance of the observed fea- 
ture we performed the following steps: 

(1) we corrected the x^ distribution by subtracting the trend mod- 
eled with a 2nd order polinomial fit. The value at Poib subtracted 
for the pedestal is 206.7. 

(2) we computed the histogram of the resulting x^ distribution be- 
tween 0.5 and 40 days [Figure |2](b)], excluding the peak region 
(29.243-31.243 days). 

(3) we fit the distribution for x^ > 15 with an exponential func- 
tion and evaluated the integral of the best-fit function beyond 206.7. 
This integral yields a number of chance occurencies due to noise of 
2.03 X 10~^^ corresponding to a significance of the detected tem- 
poral feature of ~7.3 standard deviations in Gaussian statistics. 

As an alternative method to evaluate the significance of the 
temporal feature we produced 2000 light curves performing a ran- 
dom distribution of the observed rates in the bin times of the origi- 
nal light curve. A periodogram between 0.5 and 100 days [51673.0 
trial periods with the resolution of P'^ /{N AT)] was calculated for 
each of them obtaining a maximum x^ value of 91.5 in the whole 
sample. This corresponds to a significance of the observed period- 
icity higher than 5.5 standard deviations. 

In order to exclude the presence of systematic features in the 
periodogram, we performed the same analysis on the light curve 
of the two nearby sources (IGR J16479-4514 and GX 340+0). We 
found no evidence for any significant feature at Porb . 

In Fig. |2] (c) we show the pulsed profile folded at Porb with 
Tcpoch = 54172.4236 MJD. There is clear evidence for a phase 
consistent with null intensity, whose centroid, evaluated by fit- 
ting the data around the dip with a Gaussian function, is at phase 
0.268 ± 0.015, corresponding to MJD (54180.5±0.5) inPoA . 

The BAT light curve at the time resolution of 1 day (Fig. [3]l 
(top panel) shows no evidence for the presence of flare episodes 
throughout the period of monitoring. The average count rate is 
(6.9 ± 0.5) X 10""" counts s^^ pixeP^ and the maximum count 
rate spread is consistent within 3 standard deviations with the av- 
erage. Figure |3] (bottom panel) shows the 54-months BAT light 
curve with a bin time of Porb, excluding bins with an exposure 
fraction less than 5%. The source is detected in most of the time 
bins (30 out of 38). Adoptin g the best fit model of the INTEGRAL 
data JLutovinov et alJuOOSh . the average count rate corresponds to 
an observed 15-50 keV flux of 2.36 x 1 0~" erg cm~^ s~ \2.3 
mCrab). Assuming a distance of 9.5 kpc ( Nespoli et ahluOOSl) , the 
average luminosity is 2.6 x 10^^ erg s^^. 



Figure 2. (a): Periodogram of Swift-BKI ( 15-50 keV) data for 
IGR 116465^507. (b): Distribution of x^ values extracted around Porb, 
in the period range between 0.5 and 40 days excluding the values between 
29.243 and 31.243 days. The continuous line is the best fit obtained with 
an exponential model applied to the tail of the distribution (x^ > 15). (c): 
Swi/f-BAT Light curve folded at a period P = 30.243 day, with 16 phase 
bins. 



4 SOFT X-RAY DATA 

The source was observed by Swift-XKi on 2009 June 13 (ObsId 
00037885001, 2.3 ks net exposure). The epoch of the observation 
corresponds to an orbital phase of ~ 0.2. The data were processed 
with standard techniques with the FTOOLS in the Heasoft pack- 
age (v.6.8). Figure|4]shows the background subtracted 0.3-10 keV 
light curve, with a time bin of 300 s. 

The statistic content of the XRT data is too low for a mean- 
ingful spectral analysis. In order to convert the observed count rate 
into 0.3-10 keV flux w e have used the best fit model obtained by 
iLutovinov et al.l ( l2005h for the XMM data. We obtain an observed 
flux in the 0.3-lOkeV energy band of 1.16 x 10~^^ erg cm~^ s"'^. 
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Figure 4. XRT light curve binned at 300 s. 
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Figure 3. Top panel: BAT liglit curve with a bin time of 1 day . Bottom 

panel: BAT light curve with a bin time of Porb ■ 



This corresponds to an X-ray luminosity Lx ^ 1-2 x 10 erg s 
in the same band. 



5 DISCUSSION AND CONCLUSIONS 

We have analysed the data collected by Swift-BKi^ during the first 
54 months of the Swift mission in the region of the supergiant 
HMXB IGR J 16465^507. The source is detected at a signifi- 
cance level of ~ 14 standard deviations with an average flux of 
~ 2.36 X 10"" erg cm"^ s"^ in the 15-50 keV energy band. 

The light curve reveals a periodicity of 30.243 ± 0.035 days 
that we interpret as the orbital period of the binary system. 

By applying Kepler's third law, the semi-major axis of the bi- 
nary system is given by a?' — P^j.^, x G(M* + Mx)/47r^, where 
M* and Mx are the masses of the supergiant and compact ob- 
ject, respectively. W e adopt Mx = 1.4 Mq and M* = 27.8 Mq 
JMartins et al]|2005l for an 09.5 I star of radius K^, = 22.1 Rq). 
This yields a ~ 125 Rq ~ 6 R*. We note that the assumption of 
a stel lar type B0.5 lb wo uld lead to an estimate of a ~ 150 Rq ~ 
5 R* JSearle et al.ll2008l M* = 47 Mq, R* = 32.2 Rq). 

The folded profile shows the presence of a dip with a count 
rate consistent with no emission, that could be intepreted as a full 
eclipse. However, the width of this dip (~ 20% of Porb) is not con- 
sistent with the duration of the eclipse expected for a circular orbit 
with 6 R* semimajor axis, assuming an edge-on inclination. To ob- 



tain such a long eclipse, the system should have an eccentricity of at 
least 0.8 coupled with a high inclination. This scenario is unlikely. 
Alternatively, this wide phase interval of minimum intensity could 
be explained with an eccentric orbit and a lower wind density at 
apastron. 

Given our knowledge of both a spin period and an or- 
bital period w e can locate the source on the Corbet diagram 
JCorbetl 119860 togeth er with the known Papin and Po rb of 
other binary systems JBildsten et"ZI 1 19971 : iLiu et alJ l2006h . We 
also plot in the diagram the position of the two SFXTs IGR 
J 18483-031 1 JLevine & Co rbet 2006; Zurita Heras & Chatv 200S ) 
and IGR Jl 1215-5 952 fswan k et al. 2001; Sidoli et al. 2007 ; 
iRomano et ai]|2009l) . IGR 116465^507 sits at the boundary of the 
wind-fed OB-HMXBs (systems with OB supergiants that underfill 
their Roche-lobes) and the locus of the Be transients (Fig 



;t al.l I2OO4I 



^The characteristics of the system dLutovinov et al.l I2004L 

I2OO5I) suggested a cla ssification of this source as a SFXT 
JNegueruela et alJ 120061) . However, the outburst history of this 
source is rather scarce. After the initial discovery by INTEGRAL 
no further outbursts have been reported. Wal ter & Zurita HerasI 
( l2007h mention three episodes of enhanced hard X-ray activity ob- 
served by INTEGRAL at the limit of the instrumental sensitivity, 
that cannot be considered as real flares. The BAT 15-50 keV light 
curve during the 54 months of monitoring indicates that this is a 
faint persistent source, with an average luminosity of 2.6 x 10'^^ erg 
s~^ (at a distance of 9.5 kpc), no evidence for flaring activity and 
a weak variability with a dynamical range lower than 10. Based on 
their X-ray variability, sgHMXB can be classified as classical and 
absorbed systems ( variability factor <20) or as SFXT (variabil- 
ity factor > 100). The timing behavior of IGR J16465-4507 sug- 
gests a wind-fed system, with a neutron star continuously accreting 
from a rather homogeneous wind. The luminosity, lower than what 
observed in classical sgHMXB (10^® — 10^^ erg s~^), can be ex- 



plained with the larger orbital separation ('- 
R* in classical systems). 



'5 R* compared to ^2 
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Figure 5. The Corbet diagram showing the neutron star Pspin vs. binary 
period Porb- Black circles are HMXBs with an OB primary, empty circles 
those with a Be one. Larger (red) symbols represent SFXTs: triangles are 
IGR J18483-0311 and IGR Jl 1215-5952, the star the newly determined 
position of IGR 116465^507. 
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